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We experimentally and theoretically investigate the collective behavior of three colloidal particles that are 
driven by a constant force along a toroidal trap. Due to hydrodynamic interactions, a characteristic limit cycle 
is observed. When we additionally apply a periodic sawtooth potential, we find a novel caterpillar-like motional 
sequence that is dominated by hydrodynamic interactions and promotes the surmounting of potential barriers by 
the particles. 

PACS numbers: 82.70.Dd, 67.40.Hf, 83.80.Hj 



Hydrodynamic interactions (HI) play an important role 
whenever two or more particles move in a viscous fluid III 01 . 
Due to their long-range nature, they govern the dynamics of 
colloidal suspensions, e.g;, during self- and collective diffu- 
sion sedimentation |@], and aggregation processes @]. 
Furthermore, HI can lead to pattern formation of rotating mo- 
tors |6] with a possible two-dimensional melting transition 0] 
and they are indispensable for the locomotion of microorgan- 
isms 00] or in the transport of fluid by beating cilia lioll . 
While in all these examples many colloids are involved, the 
effect of HI in few-particle systems has been investigated only 
recently. It has been demonstrated that HI mediate the corre- 
lated motion of a pair of colloids trapped in optical tweezers 
HHl and that they give rise to interesting collective behavior, 
e.g., periodic or almost periodic motions in time lll2ll or even 
transient chaotic dynamics in sedimenting three-particle clus- 
ters m . 

In this Letter, we experimentally and theoretically demon- 
strate how HI lead to a novel motional behavior of a colloidal 
system comprised of at most three particles. In contrast to 
the aforementioned examples, where the colloids exhibit ei- 
ther deterministic drift or Brownian diffusion, in the follow- 
ing we concentrate on a non-equilibrium system where both 
deterministic and stochastic motions are of importance. This 
work is partially motivated by a recent theoretical analysis of 
particles driven by a constant tangential force along a toroidal 
trap Q . Owing to HI, the particles first go through a tran- 
sient regime and then enter a characteristic limit cycle. Here, 
we present the first experimental confirmation of these find- 
ings. Our main objective, however, is to investigate experi- 
mentally and theoretically how the collective motion of inter- 
acting particles changes when a sawtooth potential is added 
to the constant driving force. Sawtooth potentials are an im- 
portant component for thermal ratchets studied, e.g., in con- 
nection with biological motors fiUl . Here, we demonstrate 
that, due to HI, two-particle clusters exhibit an unexpected 
caterpillar-like motion which facilitates the surmounting of 
potential barriers. This motional sequence is largely domi- 
nated by hydrodynamic interactions in the system. 

Tangential driving forces were exerted on colloidal particles 
with a single three-dimensional laser tweezer that scans a cir- 
cle inside our sample cell with the help of computer-controlled 



mirrors at a frequency fy. In contrast to high scanning speeds 
(/t > 200 Hz), where the particle motion is entirely diffu- 
sive fi^l . it has been experimentally and theoretically demon- 
strated that, at intermediate scanning speeds where the particle 
cannot directly follow the trap, a small periodic force is trans- 
ferred from the passing optical trap onto the particle fl7ll . As 
a result, the particle moves with a mean velocity vn « fj 1 ■ 
For the experimental parameters in our setup (A. = 532 nm, 
P « 200 mW, fj = 76 Hz), this yields vo w 7 /jm/s for a sin- 
gle silica sphere with radius a = 1.5 immersed in ethanol 
solution. Since the particle displacement by a single kick from 
the optical trap is estimated to be only about 0.08a, the parti- 
cle motion is rather smooth. In addition, the focus size slightly 
changes along the circle which in total leads to a variation of 
Vo smaller than 20%. This means that we can also view vo as 
the result of a constant driving force Icq = 67ir|flvo. To avoid 
wall effects which further complicate the theoretical treatment 
of HI, the focus of the laser beam was about 40 fjm above the 
bottom plate of our sample cell. The particles were illumi- 
nated with a white light source and imaged to a CCD camera, 
which was connected to a PC where images were compressed 
and stored. Particle trajectories were obtained offline by a 
particle tracking algorithm To enhance the electrostatic 
coupling between the charged spheres, no salt was added to 
the suspension. 

A sequence of typical snapshots (At = 4.8 s) of three par- 
ticles in the periodic limit cycle driven counter-clockwise on 
a circular ring with 9.86 radius is shown in Fig.^ The 
two-particle cluster at the top [Fig. E a )] catches up with a 
preceding third particle [Fig.^b)] until they form a triplet for 
a short time. Since the middle particle in this triplet is most 
efficiently screened from the fluid flow, it pushes the front- 
most particle ahead so that the two front particles leave the 
last one behind [Fig.^c)]. Due to drag reduction by drafting, 
this cluster is then again catching up with the single particle 
[Fig-Ed)]. We also investigated the transition from the unsta- 
ble particle configuration, where the particles had originally 
the same distance, into the periodic limit cycle and found very 
goo d agreement with theoretical predictions (data not shown) 
11411 . It should be also mentioned that the limit cycle described 
above is entirely absent when repeating the experiment close 
to a surface. We believe that this is a result of the rigid bound- 
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FIG. 1 : Snapshots of a video sequence flUl describing the charac- 
teristic limit cycle of three colloidal particles (bright) driven along 
a toroidal trap in counter-clockwise direction. The time difference 
between the pictures is 4.8 s each. The optical trap is blocked with 
optical filters. 



ary that alters HI close to surfaces |2f, 

After having demonstrated hydrodynamic effects in the 
presence of a constant driving force, we now want to study the 
cooperative particle motion in the presence of a more compli- 
cated force profile. In addition to the constant driving force 
ko > 0, we apply a sawtooth potential V s (x) with period L, 
where x is the arc-length coordinate along the circumference 
of the trap. In the first segment of length L\ , V s (x) exerts an ad- 
ditional force k s \ > on the particle, thus enhancing the drift 
motion, whereas in the second segment of length L,2=L — L\, 
the force — k S 2 < counteracts ko. As a result, the particle is 
moving in a tilted sawtooth potential V(x), as illustrated in the 
inset of Fig.|2ja). The single-particle motion is described by 
the one-dimensional Smoluchowski equation from which the 
probability distribution p(x,t) for finding the particle at posi- 
tion x in the tilted sawtooth potential is calculated. Since the 
particle moves in a circle, p(x,t) is stationary. After solving 
the Smoluchowski equation on the two segments using appro- 
priate matching conditions at x = 0, L and x = L\ , we calculate 
the constant probability current j and finally arrive at the mean 
particle velocity (v) = jL lEUl : 
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Here, vo = ko / (6nr[a) is the particle velocity without sawtooth 
potential, and h = L-ijL describes the asymmetry of the poten- 
tial. C = koL/ (ksT) is the energy dissipated on the length L 
relative to the thermal energy k%T. It is linked to the conven- 
tional Peclet number Pe = Ca/L. Since C ^> 1 in our case, 
one might naively expect a purely deterministic motion. In 
general, however, the type of motion depends on the value 
q = ksi/ko which serves as a measure for the amplitude of the 
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FIG. 2: (a) Normalized single-particle velocities as obtained from 

experiments (□) and from the analytical expression Q ( ). Inset: 

Schematic representation of the tilted sawtooth potential, (b) Cor- 
responding two-particle-cluster velocities from experiments (□) and 
numerical simulations including HI (A) and without HI (•). For 
comparison, the dashed curve of (a) is replotted. 



sawtooth potential. For q < 1, the motion is indeed purely 
deterministic and determined by the first term in the second 
line of Eq. At q = 1, the net force acting on the particle 
in segment L? vanishes (ko — k<a = 0), and the colloid moves 
entirely stochastically until it resumes the drift motion in seg- 
ment L\ . Atq> 1, it even experiences a potential barrier that 
gives rise to a stick-slip motion. Fig.|2J a ) (dashed line) shows 
the result of Eq. Q as a function of q for the set of experimen- 
tal parameters specified below. Note that, in the limit C — > °°, 
the Boltzmann tail at q > 1 vanishes and (v) becomes zero at 
q = l. 

To realize the situation described above experimentally, we 
weakly modulated (< ±12%) the intensity of the scanning 
optical tweezer. This was achieved with an electro-optical- 
modulator controlled by a function generator that was syn- 
chronized with the scanning motion of the laser focus. Be- 
fore discussing hydrodynamic coupling of particles in such a 
situation, let us briefly demonstrate that our experimental ap- 
proach allows to simultaneously apply a constant drift force 
and a quasi-static periodic potential to the colloids. The up- 
per inset of Fig.|3]shows the position dependent velocity v(x) 
of a single particle (averaged over 200 ms each) determined 
from its trajectory. From this, we calculate the position de- 
pendent force k(x) = 67tr|av(jc) and obtain the energy dissi- 
pated by the particle, Wdi ss (x) = — [q k(x')dx' (Fig. [3}, which 
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FIG. 3: Experimentally determined dissipated energy of a particle 
in a sawtooth potential with a constant driving force. Upper inset: 
position-dependent velocity of a single particle. Lower inset: saw- 
tooth potential as obtained from the dissipated energy (□) and as cal- 
culated from the intensity variation along the contour of the toroidal 
trap (---). 

is essentially identical to the tilted sawtooth potential V(x) 
ll22ll . After subtraction of the linear contribution kox corre- 
sponding to the mean particle velocity vo measured without 
intensity modulation, we finally obtain the spatially periodic 
potential acting on the particles. The open symbols of the 
lower inset in Fig.|5]show the corresponding sawtooth poten- 
tial. On the other hand, taking into account the size of the 
laser focus and that of the particles (both effects lead to some 
rounding of the edges in the potential), we can calculate from 
the sawtooth-shaped modulation of the laser intensity the ef- 
fective potential acting on the particles lElll . It is shown as 
dashed line and demonstrates the good agreement with our 
measurements. The amplitude of the potential corresponds 
to an intensity modulation amplitude of AP « 20 mW that is 
consistent with the measured intensity variation. 

First, we experimentally determined the normalized mean 
particle velocity (v) /vo of a single particle [see Eq. Q] as a 
function of the amplitude of the sawtooth potential, the latter 
being proportional to q (with the proportionality constant as 
fitting parameter). The results are plotted in Fig. 0a) as open 
symbols and were obtained for the following experimental pa- 
rameters: vo = 7.24 /jm/s and 8 = L2/L m 2.7 /jm/10.3 /jm 
which yields k = 245 • 1(T 15 N and C ps 610. The dashed 
curve in Fig. a) is the result of a least-mean-square fit to 
Eq. that shows excellent agreement with the data when L2 
as a free fitting parameter assumes the value L2 = 1 /jm and 
thus 8 = 0.1. We attribute the deviation from the experimental 
L2 to the difference between the experimentally realized and 
the perfect sawtooth potential as mentioned above. 

Next, we investigated the case where three particles were 
driven along the sawtooth potential. Similar to the case 
V s (x) = 0, the particles change their relative distances as a 
function of time (for video sequences see Jlji). In contrast 
to Fig. [2 however, we do not observe the same periodic limit 




FIG. 4: Sketch of the motional sequence of a cluster comprised of 
two particles. Due to HI and electrostatic interactions, a caterpillar- 
like motion is observed which facilitates the surmounting of potential 
barriers. 



cycle but a novel type of collective motion induced by the 
sawtooth landscape. Figure 0b) shows the average velocity 
of two-particle clusters in a tilted sawtooth potential as deter- 
mined from experiments (open symbols) and from Brownian- 
dynamics simulations (solid triangles) with the above param- 
eters. In our image analysis, we defined a cluster as a pair 
of particles with center-to-center distance < 5a, correspond- 
ing to arc-length distance < 0.7L. Obviously, the two-particle 
velocity is larger than the corresponding one-particle veloc- 
ity which is again a result of the reduced hydrodynamic fric- 
tion. What is actually surprising, is the fact that even when 
the single-particle velocity drops to zero around q > 1 (be- 
cause the particle becomes trapped in a potential well), the 
cluster velocity varies only by about 20% in that g-range. This 
clearly demonstrates that the surmounting of potential barriers 
is largely facilitated for particle clusters compared to a sin- 
gle colloid. The reason for this behavior is due to HI that 
lead to an interesting, caterpillar-like motion of the clusters as 
described in the following. Assume particle "2" [Fig. 0a)] 
drifts into a potential well which is already occupied by par- 
ticle "1". Due to a combined effect of electrostatic repulsion 
and HI, particle "2" will push particle "1" over the barrier 
[Fig. 0b)]. This in turn causes a hydrodynamic drag that 
particle "1" exerts onto particle "2" [(Fig. |^c)] which pulls 
particle "2" across the barrier [Fig.0d)]. The motional se- 
quence just described is reminiscent to that of a caterpillar 
which first stretches out, adheres at the front, and then pulls 
the tail towards the head. Once such a mode is initiated (e.g., 
by a thermally induced process), the outlined motion may last 
over several periods L until the particles become trapped in 
separate potential minima due to thermal fluctuations of their 
distance. Then, the motion stops until the situation described 
in Fig. 0a) is initiated by a spontaneous hop of particle "2" 
into the potential well occupied by particle "1". 

This motional pattern of two-particle clusters was also 
observed in our Brownian-dynamics simulations. The 
electrostatic inter-particle repulsion is described by a 
screened Coulomb potential Q V lep (r) = [Ze Ka /(l + 
-Ka)] 2 X^k^Te^ Kr / r with Debye length k -1 ps 300 nm, where 
the effective particle charge is adjusted to Z ps 8000 (taking 
A-b ~ 0.7 nm for the Bjerrum length in water) so that the 
center-to-center distance between the particles does not go be- 
low about 3a as observed in the experiment. The Brownian- 
dynamics simulations including HI were carried out using a 
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predictor-corrector-type integration scheme which is first or- 
der in the time step Af I2H1 . For particle separations r > 3a, 
the mobilities describing the HI and thus the mutual coupling 
of all particles can be well approximated by the Rotne-Prager 
tensor |2], which is the far- field expansion up to order 1/r 3 
(for details, see Ref. El ). 

To demonstrate the crucial role of HI for the enhanced clus- 
ter motion, we performed numerical simulations without HI. 
The mean particle velocity, defined as the average over all 
particles and times, coincides with the single-particle curve 
in Fig.|2]when plotted as a function of q. This is clear since, 
due to the average over all particles, the electrostatic forces 
cancel each other (actio = reactio) and thus, on average, the 
particles move independently. The mean velocities of the par- 
ticle clusters, however, deviate from the single-particle veloc- 
ities, as illustrated in Fig. |2jb) (solid circles). Nevertheless, 
as expected, the average cluster velocities are always below 
the corresponding curve where HI are included. For q > 1, 
the two-particle-cluster motion without HI consists of singular 
events where particle "2" moves down the seceding flank of 
the tilted sawtooth potential and pushes particle "1" over the 
potential barrier [Fig.|4la,b)]. In contrast to the caterpillar-like 
behavior, the collective motion stops here since the drag due 
to HI is missing. Then, another thermal activation is necessary 
to induce the sequence again fl9ll . 

Finally, we point out that the short-range electrostatic re- 
pulsion between the particles is only of minor importance for 
the collective, caterpillar-like motion as described above. It 
is rather the long-ranged HI that dominate the force which 
pushes the front particle in Fig. |4jb) over the barrier (note 
that HI decay asymptotically as 1/r whereas the domi- 
nant term in the repulsive electrostatic force is exponentially 
damped, F rep (r) = -V,V rep (r) °c e~ r / x ° jr with X D < a). The 
caterpillar-like motion was also obtained in simulations where 
the strength of the electrostatic interaction potential was re- 
duced by a factor of 50 so that the minimum distance to con- 
tact of the particles was 0.1a lliUl . 

The present study demonstrates that HI strongly dominate 
the motional behavior of driven colloidal particles. In the 
presence of a constant driving force, we experimentally con- 
firmed the recently predicted characteristic limit cycle. If in 
addition a static periodic potential is applied along the toroidal 
trap, we find that colloidal clusters perform a caterpillar-like 
motion which facilitates the surmounting of large potential 
barriers. This novel type of motion which is predominantly 
triggered by HI is also confirmed by numerical simulations. 
In addition to a better microscopic understanding of HI in 
few-particle systems, our results are also of relevance for the 
motion of interacting particles in thermal ratchets, which are 
considered as models for biological motors. Indeed, recent in- 
vivo experiments demonstrate that the speed of coupled motor 
prot eins is increased compared to the speed of a single motor 
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